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I. INTRODUCTION 

Spin chains are gaining increasing attention as natu- 
ral candidates of quantum channels 0, 0, 0, 0, ■ An 
unknown quantum state can be prepared on one end of 
the chain and then transferred to the other end by sim- 
ply employing the ability of the chain to propagate the 
state by means of its dynamical evolution. This proce- 
dure does not require any gating and therefore can be 
implemented without any need for a modulation of the 
couplings between the spins. Especially this last aspect 
may be of importance for solid state quantum informa- 
tion. Several aspects of quantum communication using 
spin chains have been already obtained starting from the 
original proposal by Bose The fidelity of the trans- 
mitted state can be increased if the sender and the re- 
ceiver can encode/decode the state using a finite number 
of spins Q. Even perfect transmission can be achieved 
if the exchange couplings of the chain are chosen appro- 
priately or if it possible to perform measurements on 
the quantum spins of the chains ,6|. Moreover, it was 
recently shown that a phase covariant doner can be real- 
ized using an unmodulated spin network Q and that the 
spins chain approach to quantum communication seems 
to be realizable in solid state devices with present day 
technology 

Another interesting aspect which is emerging recently 
is the interconnection between quantum information and 
condensed matter theory Examples are the study of 
non-local (quantum) correlations of spin systems in a va- 
riety of situations (see |lfjj | and references therein) or re- 
formulation of the density matrix renormalization g roup 
in the framework of quantum information theory 
In this paper we would like to further explore the use 
of concepts born in quantum information for the char- 
acterization of spin chains. Here we propose to use the 
information capacities of quantum channels as a tool to 
characterize some aspects of dynamical correlations in a 
spin chain. By calculating the capacities of the channel, 
obtained by identifying two separate sections of the chain 
as the extremes of a communication line, one can in fact 
get some information about the strength of the corre- 
lations among the interconnecting spins. In some sense 
this is equivalent to assigning to any couple of subsets of 
the spins in the sample the values of the corresponding 
capacities in order to create a "road map" of the infor- 



mation fluxes in the system. 

Our results apply for the whole class of Hamiltonians 
for which the total magnetization along a fixed direc- 
tion is conserved. In the present work we confine ourself 
to the case in which one bit is transferred through the 
chain. This situation corresponds to study the sectors in 
which only few spins are up. We believe, however that 
the present approach can be further extended to other 
spin sectors. 

The paper is organized in two distinct parts. In the 
first one we introduce the spin chain communication lines 
and we discuss some models in details (Sec. [HJl. I n tne 
second part (Sec. IHIfl instead we give a brief overview of 
channel capacities and, by studying the qubit amplitude 
damping channel |T^ , we calculate the quantum capacity 
Q, the entanglement assisted capacity Ce and the clas- 
sical capacity G\ of the spin chain models presented in 
Sec. [HJ As explained in detail throughout the paper the 
above quantities are computed in communication scenar- 
ios where many "copies" of the spin system are available 
and "not" for those scenarios where the communicating 
parties make multiple uses of the "same" spin chain. The 
paper ends in Sec. IIVI with the conclusions. 



II. THE MODEL 

Given a collection of N spins coupled by means of 
a time-independent Hamiltonian it is possible to define 
a quantum channel by identifying two set of spins of 
the sample (say, sets A and B) as two quantum reg- 
isters. A first party (the sender of the message) en- 
codes some information on A and a second party (the 
receiver) tries to recover such information from B some 
later time t 

[JH0. Formally, one assumes that the 
encoding procedure takes place by initially decoupling 
the spins A from the remaining spins without disturbing 
them from their initial fiduciary state oo, preparing A 
on some input message pa, and finally allowing pa and 
Co to interact for a given time t through the Hamilto- 
nian of the system. At this point we consider the state 
PB(f) of the spins B, obtained from the spin chain state 
R(t) — U(t)(pA <8> ao)W(t) by tracing away all the the 
degrees of freedom of the system but those relative to B. 
[Here U(t) is the unitary evolution of the system]. The 
resulting Completely Positive, Trace preserving (CPT) 
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mapping 

p A -> = p B (t) = Tv {B \u{t){ PA <g> a )U\t)] (1) 

where Tr' 5 -* means trace over all the spins but B, defines 
the quantum channel we are interested in. A caveat is in 
order. Equation JQ) does not provide a proper description 
of any realistic scenario where the communicating parties 
keep on operating on their quantum registers over an ex- 
tended period of time. In fact, after a first "reading out" 
of the signal from the spin chain at time t, the state <To 
of the spins not belonging to A will change, thus mak- 
ing Eq. (JTJ unable to describe the state of B for later 
times. As a matter of fact, in our system any repeated 
(in time) manipulation of the quantum registers will nec- 
essarily introduce memory effects in the communication 
for which a satisfying quantum information theory is still 
missing (some preliminary results on quantum channels 
with memory can be found in Ref. |T^1- 

Having the previous observations in mind, Eq. Q can 
still be used to provide a quite complete characterization 
of the correlations between the registers A and B. Tech- 
nically the map M. is defined in the Hilbert space Ha of 
the configurations of A and depends on t, Hi^t. an and on 
the choice of A and B. A common approach 0,1213 used 
to quantify the correlations between the sets A and B of 
the chain is then to consider the fidelity associated with 
M. , by computing the average fidelity between the input 
states pa and their output counterparts ps(i). Here we 
observe that a more detailed characterization of such cor- 
relations can be obtained by treating M. of Eq. as the 
CPT map of a real memoryless channel. The correlations 
between A and B can then be analyzed by means of the 
information capacities [12, Uj, UH Ha, LLll UM associated 
with M.. These quantities are related to the minimum 
"amount of redundancy" over multiple uses of the chan- 
nel M. needed to achieve perfect message transmission, 
i.e. unitary fidelity: higher values of the capacities corre- 
spond to stronger correlations between input and output 
states. Depending on the nature (e.g. classical or quan- 
tum) of the information propagating through the channel 
one can define different capacities of M. and each of them 
are obtained by maximizing over all possible coding and 
decoding strategies which act over multiple channel uses 
(see Sec. Illlf) . From the discussion that follows Eq. 
it should be evident that the capacities of the map M 
might not provide a proper description of the communi- 
cation performances of the spin chain. As a matter of 
fact, they only account for those scenarios where many 
parallel "copies" of the same system are simultaneously 
operated by the two communicating parties but fail to de- 
scribe those scenarios where the communicating parties 
make successive multiple uses of the "same" spin chain. 

A. Solvable models 

Examples of a spin chain channel where the informa- 
tion capacities of M. can be solved exactly are provided 



by a chain of 1/2-spins coupled through a ferromagnetic 
Heisenberg interaction . The system Hamiltonian is 

N 

H = -Y,*Jid (44 + 44 + ^44) - E nB *4 ,(2) 

where the first summation is performed on the nearest- 
neighbor spins of the chain, a l x y z are the Pauli operators 
associated with the ith spin, jy are coupling constants, 
7 is an anisotropy parameter, and -B, are associated with 
externally applied magnetic fields. 

In the following we will assume the input set A and 
output set B of Eq. to contain, respectively, the first 
k and the last k spins of the chain, and that all the spins 
but A are initially prepared in the same eigenstate | J, ) of 
a z . [Here | J, ) and | f ) are the eigenstates of a z associ- 
ated respectively with the eigenvalues —1 and +1]. If the 
two communicating parties are allowed to use the entire 
Hilbert space of their quantum memory, Eq. Q yields 
a A:-qubit channel which, in general, is too complex to 
be analyzed on a complete basis. In the present paper 
we consider thus a simplification of this scenario where 
the sender and the receiver use their k spins to encode 
a single logical qubit. This approach is clearly inefficient 
from an informational theoretical point of view, but on 
the positive side, the resulting maps can be treated ana- 
lytically. 



1. Encoding one logical qubit with one-spin up vectors 
As in Refs. 0, we introduce the one-spin up vector 

ii) = iu--m---i>, (3) 

which for j = 1, • • • , N represents the state of the chain 
where the jth spin is prepared in the eigenstate | | ) and 
the other N — 1 ones in | J. ) . Suppose that at time t = 
the sender prepares her/his spins in 

\^) A = a\^) A +mi)A, (4) 

where a, (3 are complex amplitudes, | -IJ- )a is the state of 
A with all spins down, and \4>\)a is a given normalized 
superposition of \j) with j referring to spins of the input 
memory A, i.e. 

k 

\^i)a = J2 c ^)- ( 5 ) 

Since the Hamiltonian of Eq. J2J commutes with the total 
spin component along the z direction, one can show that 
at time t the whole chain is described by 2] 

JV k 

\^(t))^a\^)+(3j2J2 c ^r4t)\f) , (6) 
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with | JJ. ) the state of the chain with all spins down and 
with 

UAt) = (j\e- lHt/h \s) (7) 

[see, for instance, Ref. Q for the explicit functional de- 
pendence of fj, s (t) from the evolution time t\. According 
to Eq. Q the state of B at time t is finally obtained from 
Eq. |JBJ by tracing over all the remaining N — k spins, i.e. 



PB (t) = (M 2 + (i-77)|/3| 2 )|4M4l + ^mM0' 



with 



N 

E 



j'=N-k+l j=l 



(8) 



(9) 



In Eq. JSJ the two orthonormal vectors | JJ. )b and \4>[}b 
are, respectively, the state of B with all spin down, and 



N 



(10) 



Apart from an irrelevant unitary transformation [T^ the 
map associated with _Ps(i) of © is a qubit amplitude 
damping channel 0, Q of efficiency rj which acts on the 
orthonormal basis {| JJ. )b, \4>'i)b}- For the sake of clarity 
let us identify | 4J- )a an d | JJ- )b with the same logical 
qubit state |0) and \cf)i)A, \(j>[) b with |1). In this notation 
we can express the input state | \E'} ^4 and the output state 
PB{t) as density matrices p and p' of the same qubit 
Hilbert space Ha- Equation JHJ becomes thus 



p' = V^p) 



(11) 



with the amplitude damping map characterized by 



the Kraus operators ^3 



Ao = |0)<0|+V*?|1><1| , 



A x = vT^|0)<l| 



(12) 



Equation (|11() describes a quantum channel in which the 
logical information of A represented by the coefficients a 
and (3 of Eq. (0J , is transferred to the output memory B 
with an accuracy which can be estimated by calculating 
the capacities of the map T> n . The calculation of the 
capacities of T> v is presented in Sec. IIII Al 



In other words, the sender still uses the state with no 
spin up to "transfer" a, but now a selected superposition 
of two-spin up states is employed to "transfer" /?. In this 
case one can show that the output state of the memory 
B is 



PB (t) = (M 2 + 772|/3| 2 )l4M4l + ml/?l 2 l02M 
+ -m)\P\ 2 °B , 



(14) 



where |0' 2 )s is a superposition of two-spins up state of 
B and erg is a density matrix of one-spin up states of B 
whose eigenvectors are orthogonal with respect to | JJ-) B 
and \4>2)b (see Appendix lAl for details). For r/i + r\% = 1 
the mapping (|14|l reduces to a qubit amplitude damping 
channel T> v with quantum efficiency t] — r\\. However, in 
the general case, Eq. (|14fl is slightly more complex. In 
fact for t}\ + 772 < 1 the component \4>2)a °f t ne input 
state undergoes to three possible processes: with proba- 
bility 7/i it is rotated into the output state |(t3 2 )b; with 
probability 7/2 it is damped to | JJ- ) b and finally with 
probability 7/3 = 1 — 7/1 — t/ 2 it is transformed in the den- 
sity matrix ob- In this respect, the map l|14(l is similar 
(but not equal) to a channel that, with probability 7/3, 
decoheres the input and transforms \4>2}a into o~b while 
with probability 1— 7/3 applies to it an amplitude damping 
channel transformation of quantum efficiency 7/1/(1 — 7/3). 

A compact description of l|14|) is obtained by identify- 
ing I JJ. )ai I JJ- )b with the logical qubit state |0) and 
\4>2)a, \4^2)b w hh |1). With this notation, the transfor- 
mation l|Ti|) can be expressed as a two-parameter CPT 
map, 



P' = T ni,m(p) 
with Kraus operators given by 

^0 = 
A 1 -- 

A 1+i = 



(15) 



m|l)(l| 

m\o)(i\, 



(16) 



Here \d) are the eigenvectors of <tb associated with the 
eigenvalues Q > 0: according to Eq. (|A8|) they are or- 
thogonal with respect to |0) and |1) and there are at 
most k of them. 

The capacity of the channel % n , m is derived in 
Sec. lHIBl 



2. Encoding one logical qubit with two-spin up vectors 

Consider the case where in Eq. the vector l^i)^ 
is replaced by a normalized superposition \(j)2)A of states 
\j, £} where the j-th and ^-th spins of the chain are in 
t ) while the remaining are in | | ), i.e. 

k 

\h) A ^ £ djMQ ■ (i3) 

j>£=l 



III. CHANNEL CAPACITIES 

The quantities we are interested in this paper are the 
quantum capacity Q, the classical capacity C and the en- 
tanglement assisted capacities Ce and Q e ■ The quantum 
capacity Q measures the maximum amount of quantum 
information that can be reliably transmitted though the 
map M. per channel use [15^ . Intuitively, this quantity is 
related with the dimension of the largest subspace of the 
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multi-uses input Hilbert space which does not decohere 
during the communication process. The value of Q (in 
qubits per channel uses) can be computed as 

Q = sup Q n /n (17) 

n 

with 

Q n = max \s(M® n (p))-S((M® n ® W)($))| (18) 
pen's" l J 

On one hand, the sup in Eq. (|17fl is evaluated over n 
parallel channel uses [2CJ, where the channel map is de- 
scribed by the super-operator which transforms the input 
states p of U® n into the output states M® n (p). On the 
other hand, for fixed values of n, the maximization in 
Eq. (fl^|l is performed on all possible input density matri- 
ces p £ 7Y® n . The quantity in the bracket is the coherent 
information of the channel [Titl^. S(p) = — Tr[<9lo g 9 p] 
is the von Neumann entropy and $ is a purification [l2( 
of p 6 H®" defined in the extended space obtained by 
adding an ancillary space Ttanc to H® n . 

The classical capacity C gives the maximum amount 
of classical information that can be reliably transmitted 
through the channel per channel use: here, in the multi- 
uses scenario, the goal is to identify the largest set of 
orthogonal input messages which remain distinguishable 
(i.e. orthonormal) during the propagation. In this case, 
the system is not required to preserve the phases of super- 
positions of different messages: C characterizes thus the 
ability of the channel M. in preserving occupation num- 
bers but not its decoherence effects on the transmitted 
signals. 

As in the case of Eq. (|17fl . in calculating the classical 
capacity it is necessary to perform a maximization over 
multiple uses of the channel Ai, i.e. [l6( 

C = sup C n /n (19) 

n 

where C n is the classical capacity of the channel which 
can be achieved if the sender is allowed to encode the 
information on codewords which are entangled only up 
to 7i-parallel channel uses. The value of C n is obtained 
by maximizing the Holevo information |22| at the output 
of n parallel channel uses, over all possible ensembles 
{£,k,Pk}, i-e. 

C n = max (s(^»°(p))-Vy(M«"( ft ))), 

(20) 

with {£fc} probabilities and p = J2k^ k P k ^ ne average 
message transmitted. The optimization lfH?|) could be 
avoided if the additivity conjecture of the Holevo infor- 
mation is true 0,|2^|: in this case in fact the optimal en- 
sembles pk\ which achieve the maximum in Eq. I|2U|) 
are separable with respect to the n parallel uses and C 
coincides with C\. 



The last capacities we consider are the entanglement- 
assisted classical capacity Ce and its quantum counter- 
part Qe [l7lll8l| . These quantities give, respectively, the 
maximum amount of classical or quantum information 
that can be sent reliably through the channel per chan- 
nel use, assuming that the sender and the receiver share 
infinite prior entanglement. To calculate Ce it is not re- 
quested to perform a regularization over parallel channel 
uses as in the case of Eqs. (fT7)l and l(T9")l. Here instead 
one has to maximize the quantum mutual information 
for the single channel use [13, EE UM , i-e. 

C E = max is{p) + S(M(p)) - S((M ® l Qnc )($)))(21) 

where now $ is a purification of the input message p G TL. 
The entanglement assisted quantum capacity can then be 
obtained as Qe — Ce/2 by means of quantum telepor- 
tation [25| and superdense coding [26|. The relevance 
of Ce relies on the fact that this quantity gives a sim- 
ple upper bound for the other capacities. Moreover it is 
conjectured to provide an equivalence class for quantum 
channels [Tsl ]. 

In the next sections we will focus on the capacities of 
the CPT maps associated with spin chain models intro- 
duced in the Sees. Ill A ll and III A 21 



A. Capacities of the amplitude damping channel 

In this section we analyze in detail the qubit amplitude 
damping channel. In particular we calculate its capacity 
Q showing that in this case the maximization l|18fl over 
parallel uses is not necessary. Moreover, we derive the 
capacities Ce and C\ (classical capacity achieved with 
unentangled codewords) which were originally given in 
Ref. 18] without an explicit derivation. 

The map T> v is completely characterized by the pa- 
rameter 77 and can be seen as an instance of the map 
£ v associated with the lossy Bosonic channel [27], |28|, [23 
(see App. [E]for details). In particular T> v has the useful 
property that by concatenating two amplitude damping 
channel with quantum efficiencies rj and 77', one obtains a 
new amplitude damping channel with efficiency 7777', i.e. 

2V(X>„(p))=ZW(p)> (22) 

which applies for any input state p. 

Quantum capacity:- An important simplification in 
calculating the quantity (|17f) derives by introducing the 
following representation of the CPT map V v , 

V v (p) = Ti c [V(p(g>\0)c(0\)V^}, (23) 

obtained by adding to the Hilbert space Ha of the in- 
put logical qubit A an auxiliary Hilbert space Tic, an d 
introducing the unitary operator V which in the compu- 
tational basis {|00), |01), 1 10) , |11)} of H A ® Tic is given 
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by the 4x4 matrix, 



V = 






Vo 















T) ^/rj 
1 



(24) 



In Eq. (|23(l . Trc[- • • ] is the partial trace over the elements 
of the auxiliary space Tic- The complementary channel 
T> n of T> n is defined by replacing this operation with the 
partial trace over 7i A , i-e. |30| 



V v {p) = Tr A [V(p®\0)c(0\)Vl] 



(25) 



Upon a swapping operation S which transforms A into 
C and vice-versa, the transformation of Eq. (|25|l can be 
seen as a mapping from the Hilbert space TLa into itself. 
Moreover, by direct calculation one can verify that 



V n { P ) =SV 1 - n {p) s . 



(26) 



Using the composition rule of Eq. 122|) it is thus possible 
to show that, for 77 ^ 0.5, one has [Mj 



T> n (p) = S £>(!-,,)/,, (£>„(/>)) S . 



(27) 



[The quantum capacity in the case rj < 0.5 is simple to 
compute and will be discussed at the end of this section] . 
This relation shows that for the channel V v is degrad- 
able, i.e. there exists a CPT map defined by the super- 
operator S I'(i_ f) )/ f) (- • • ) S which connects the output 
state T> v (p) with the output state T> v ( p). According to a 
theorem proved by Devetak and Shor |30[ , this condition 
guaranties that the sup in Eq. I|17fl is achieved for n = 1 
(single channel use) : in other words, Eq. (|27|l guarantees 
the additivity of the coherent information for the channel 
T> ri . For r\ ^ 0.5 the quantum capacity of derives thus 
by solving the maximization l|18|) for n = 1. 

In the computational basis {|0), |1)} of Ha, the most 
general input state of the map T> v can be parametrized 
as follows 



l-p 7* 
7 P 



(28) 



where p G [0, 1] is the population associated with the 
state |1) and | — / 1 y(X — p)p is a coherence term. A pu- 
rification $ = |$) ($| of p is then obtained by introducing 
an ancillary qubit system TL a nc and considering the state 



\<S>) = ^T^ p|0)<8|J2o>+VP|l>®|fli>, (29) 
with \Rq 1) unit vectors of H an c such that 



(Ro\Ri) = 7 A/(1-P)P 



(30) 



From the Kraus decomposition (|12() one can verify that 
the map T> v transforms this state into the output 



v v (p) = ( 1 - riP ^ 



(31) 



The matrix (|31(l has eigenvalues 

A± (v) = (l ± v/(l - 2 77p)2 +4t7|7| 2 ) /2 , (32) 
which gives an output von Neumann entropy equal to 



S(V v (p)) = H 2 (\ + (r,)) 



(33) 



with H2 the binary entropy function defined in Eq. 1C4|) . 
Analogously, by applying the map {V v ® l a nc) to $ we 
get the state described by the 4x4 complex matrix of 
Eq. IjClfl . which has entropy equal to 

S((V V ® lanc)(*)) = H 2 (X+(l-v)) ■ (34) 

The quantity we need to maximize to obtain Q is hence 
given by 

J(p, | 7 | 2 ) = H 2 (\+(r))) - ff 2 (A+(l - V )) (35) 

which depends from p only through the parameters p £ 
[0, 1] and |7| 2 e [0, (1 - p)p). As discussed in App.EJthe 
maximization (|17l) is achieved by choosing 7 = 0, which 
gives a quantum capacity equal to 

Q= max { H 2 (r ] p)-H 2 ((l-7 1 )p)\ . (36) 

For any given 77 ^ 0.5 this expression has been solved nu- 
merically and the results are plotted in Fig. ^ In Fig. |3 
instead, we have reported, as function of the parameter 
•q, the optimal value of the population p which provides 
the maximum of the right-hand-side term of Eq. I|36|) . 

Let us now consider the low transmissivity regime 
•q < 0.5. In this case, the same non-cloning argument 
given in Ref. j3^ for the erasure channel and in Ref. 
for the lossy Bosonic channel, can be used to prove that 
the quantum capacity Q of T> v nullifies. An alternative 
proof of this fact, can be obtained by noticing that the 
composition rule l(2*2l implies that the quantum capacity 
of the channel T> v is an increasing function of the trans- 
missivity r). The thesis then follows from the fact that 
for j] = 0.5 the right-hand-side term of Eq. (|36|l nullifies. 

Entanglement assisted capacity:- To calculate this 
capacity we need to perform the maximization of Eq. I|21|) 
for all possible input states p. The quantum mutual in- 
formation of the channel can be obtained by summing the 
coherent information J(p, \"f\ 2 ) of Eq. l-iol) to the input 
entropy of the message, i.e. using the parametrization 
introduced in the previous section, 

I(p,\j\ 2 )^ J( P ,| 7 | 2 ) (37) 
A + V(l-2p)2+4|7| 2 \ ^ 



According to the property 1) of App. IC li the last term 
in the right-hand-side of this expression is a decreasing 
function of I7I 2 , i.e. it is maximum for 7 = 0. From the 
previous section, we know that the same property applies 
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FIG. 1: Plot of the capacities (in bits per channel uses) of 
the qubit amplitude damping channel T> v as a function of the 
noise parameter n. The curves Q, Ce and C\ represent respec- 
tively, the quantum capacity, entanglement assisted classical 
capacity, and the classical capacity achievable with unentan- 
gled encodings: they have been obtained by solving numeri- 
cally the maximizations of Eqs. 1361 . 1381 and 14211 . For 77 = 1 
(no noise limit) both Q and Ci give one qubit per bit for chan- 
nel use, while as a consequence of the superdense coding effect 
prjl Ce gives two bits per channel uses. For 77 — 0.5 we have, 
respectively, Q = 0, Ce = 1, and Ci = 0.4717 (this result is 
in agreement with what found in Ref. H^|). As a consequence 
of the non cloning theorem 27, 32] for r/ < 0.5 the quantum 
capacity Q nullifies. The curve Qe represents the entangle- 
ment assisted quantum capacity and is obtained by simply 
dividing by two the values of Ce- The classical capacity C of 
the channel is lower bounded by Ci and upper bounded by 
1 (maximum entropy of a bit) and by Ce- If the additivity 
conjecture (23I applies to this channel then C = C\. 



also to J(p, I7I 2 ): we can thus conclude that, for any 
p G [0,1], the function I(p,\ r y\ 2 ) achieves its maximum 
value for 7 = 0. In other words, we can compute the 
entanglement assisted capacity Ce as 

C E = max { H 2 (p) + H 2 ( V p) - H 2 ((l -r))p)}- 

(38) 

This maximization can now be solved numerically: the 
resulting plot is given in Fig. ^ The optimal p's that 
saturate the maximization of Eq. I|38|l are reported in 
Fig.H 

Classical capacity with unentangled encodings:- In 
the case of the lossy Bosonic channel £ v with constrained 
input average photon number the additivity property of 
the Holevo information has been proved 28]. Unfortu- 
nately, this derivation relies on some specific properties 
of the coherent input state of the Bosonic channel: even 
though £ v and T> v are strongly related it is hence difficult 
to use the result of |28| to establish the additivity conjec- 
ture for the qubit amplitude damping channel (see also 
App. [BJ). Here we will not discuss further this problem 
and simply we will focus on the capacity Ci, which mea- 
sures the maximum amount of classical information that 



FIG. 2: Plot of the optimal populations p associated with 
the state |1) which provide the capacities of the amplitude 
damping channel T> v , as a function of the transmissivity 77. 
The upper curve refers to the entanglement assisted capacities 
Ce and Qe'- for given n it has been obtained by finding that 
value of p which maximizes the function at the right-hand-side 
of Eq. (0HJ. The lower curve refers to the classical capacity 
Ci through Eq. (1421 . Finally the intermediate curve refers 
to the quantum capacity Q through Eq. I|36|l : this function is 
not defined for n < 0.5 because for this value the Q is null (see 
Fig. 0. For n ^ 0.5 the optimal population p of the quantum 
capacity is bigger than the corresponding population of C\ 
and lower than that of Ce- 



can be reliably transmitted using only encod ing s that are 
not entangled over parallel channel uses 0, 0, . The 
quantity C\ is a lower bound for C and coincides with 
it provided the additivity conjecture holds. The classical 
capacity C\ can be calculated by solving the maximiza- 
tion of Eq. (|20|) in the case of n = 1. Consider the ensem- 
ble of messages where with probability the channel is 
prepared in the input state 



Pk 



l -Pk it 
Ik Pk 



(39) 



with pk and 7fc defined as in Eq. Q28[) . Using the result of 
the previous section we can express the associated Holevo 
information as 



1 + y/(l- 2 V p)2+4r 1 bf 



(40) 



1 + v /(l-2ryp fc ) 2 +4?7|7 fc | : 



where now p = J2k £kPk and 7 = J2k ^lk are the pa- 
rameters (|28|l associated with the average input message 
p = J2k^ k P k - According to Eq. l(2Tj|) the capacity C\ 
is obtained by maximizing x over all possible choices of 
Pk, 7/c and To solve this problem we first derive an 
upper bound for C\ and then we show that there exist 
an encoding p^, jk and which achieves such an upper 
bound. 
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From the property 1) of the binary entropy (|C4fl given 
in App. IC II we can maximize the first term in the right- 
hand-side of Eq. 1)40(1 by choosing 7 = 0. Moreover, one 
has 



E(kBt i+yTT^W+w (41) 



1 + ^1-477(1-^ 



> H2 



l + v /l-4r/(l-r / )(E fe aPfc) 2 



where the first inequality derives from the property 1) of 
H2{z) and from the fact that |7fc| 2 ^ (1 — Pk)pk, while 
the second inequality is consequence of the property 2). 
Replacing the above relation in Eq. I|40(l we obtain an 
upper bound for \ which does not depend on 7^ and 
which depends on ^ and pk only through p — J^k^kPk- 
By maximizing this expression over all possible choices 
of the variable p we get the following upper bound of C\ , 



C\ max \H 2 (rjp) 
pe[o,i] I 



(42) 



-H, 



1 + y/l- 4^(1-77) p 2 



The right-hand-side term of this inequality is indeed the 
value of C\ . This can be show by noticing that for any 
p € [0, 1] and d > 1, the parameters 



= l/d p k = p 
P)P , 



Ik = e 2 ^ d y/(T 



(43) 



with k — 1, ••• ,d, produce a Holevo information \ of 
Eq. I|40|) which is coincident with the quantity in the 
brackets on the right-hand-side of Eq. The quan- 

tities (|43[) provide hence optimal encoding strategies for 
C\. On one hand, any ensemble element pk of this en- 
coding has maximum absolute value of the coherence 
term jy- this minimizes the negative term of the Holevo 
information. On the other hand, the average message 
p = ^2k^kPk has minimum value I7I, i.e. 



7 = X>™ = \/RpE eM/d / (l = . ( 44 ) 
fc=i fc=i 

which maximizes the positive contribution to the Holevo 
information. This property of the channel T> v is a 
common feature of many other channels whose classi- 
cal capacity C has been solved [2^, 0. The value of 
C\ obtained by maximizing the right-hand-side term of 
Eq. (|42|l has been plotted in Fig. ^ while the optimal p's 
are plotted in Fig. [21 



T 



B. Capacities of the channel %j liV2 

In this section we analyze in details the CPT map 
of Eq. i|15|) associated with the spin chain com- 
munication line of Sec. Ill A 21 In this case we calculate 
the capacities Q, Ce and we provide a lower bound for 
C. 

The map T Vl iV2 is described by the positive parameters 
771 and 772 of Eqs. (|A4|) and l|A5() that satisfy the relation 
Vi + V2 ^ 1- In particular, for 772 = 1 — 771, %, ltV2 reduces 
to an amplitude damping map of transmissivity 771, i.e. 



(45) 



for any input state p. Moreover the following composition 
rule applies, 



(46) 



with 77" = 771 77^ and 77 2 ' = 772 + ViV^- An interesting way 
to express T Vl )J)2 is finally provided by concatenating two 
amplitude damping channels with a CPT map V that 
transforms one of the eigenvectors of the matrix ub (say 
I Ci ) ) m to ob itself, i.e. 



(47) 



where the first amplitude damping channel T>i^ m acts 
as usual on |0) and |1), while second one V is 
instead defined on the subspace generated by |Ci) and |1). 
The main consequence of Eq. (|47l) is that any capacity of 
%ii ,V2 cannot be greater than the corresponding capacity 
of T>i- V2 . In fact, by applying the CPT transformations 
V and D 1 // 1 _ Jj2 j to the output of an amplitude damping 
channel of transmissivity 1 — 772 one can simulate the 
corresponding output of T m )tj2 . 

Quantum Capacity:- As in the case of T> v we can 
prove that T Vl ^ 2 is degradable when its quantum capac- 
ity is not null. In fact, as in Eq. 112MI) define the unitary 
operator V of the extended Hilbert space Ha ® He such 
that for any p of Ha 



T ril , l2 (p)^Tr c [V(p®\0}c(0\)V^} 



(48) 



An example of V can be obtained by introducing the 
following vector of Ha <8> He 



\$„) ee Va\Ci)A <8 ic 



i C 



(49) 



where \C%)a € Ha are the eigenvectors of ob introduced 
in Eq. I(16J) while \C%)c is an orthonormal set of states of 
He that are orthogonal to \0)c and \l)c- The state |$ CT ) 
is a purification of ob on Ha®Hc- The unitary operator 

V can now be chosen to be the identity everywhere but on 
the subspace of Ha <8> He generated by the orthonormal 
vectors {|00), |01), |10), |<J> CT )}. On this subset we define 

V to have the matrix representation 



1 i+,/m 



V 




(50) 
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with 7/3 = 1 — 771 — 772 . The complementary map T m , V2 
is finally obtained by substituting in l|48[l the trace over 
C with the trace on A From Eqs. (03 and {2| one 
can easily verify that for 771 ^ 772 the following relation 
applies for any input p, 



%ii,V2 (P) = S /D n 2 /r il ('Tm,m{p)) s 



(51) 



with S the swapping operation which transforms A in 
C. Since T> r/2 / r]1 is CPT the above equation shows that 
for r]i ^ r]2 the map % 11 , V2 is degradable: the quantum 
capacity of this channel can be hence computed from 
Eq. (|18|l for n = 1. A straightforward generalization of 
the qubit amplitude damping channel analysis shows that 
Eq. H35fl still applies by replacing A± (77) of Eq. 13211 with 



A± fai, 772) 



1 - V3P 



(52) 



1± 



2771P 
1 - 773J? 



4m H 2 
(1 - r] 3 p) 2 



The quantum capacity of T ni tV2 becomes hence 




Q 



H 



tax { (1 - 773P) 



(53) 



1 



mp 



-Ho 



1 - (1 -77i )p N 

1 - V3P 



)]} 



Notice that as in the case of Eq. i|36|) the maximization 
over the input parameter 7 of Eq. (|28|l has been satu- 
rated by setting 7 = (the proof goes as in Eq. J2f)J|). 
A plot of Q as a function of 771 and 772 is reported in 
Fig-Hby solving numerically the maximization on p. The 
above results do not apply for 77! ^ 772: in this case in 
fact T> V2 / ni is not CPT. However, a non-cloning argument 
can be used to prove that in this parameters region, the 
quantum capacity of T Vl jJ)2 nullifies |35j . 

Entanglement assisted capacity:- The analysis of C'e 
proceeds as in the case of the amplitude damping chan- 
nel. Here Eq. H38f) is replaced by 



max 



Ho 



[01] ( + ^ ~~ 

- m)p\ „ /1 - (1 - 77i)jj n 



(54) 



mp 



-Ho 



V3P 



)]} 



A numerical plot of this expression is given in Fig. 

Lower bound for C\ :- The analysis of C\ for the chan- 
nel T Vl >V2 is slightly more difficult than that for T> v as the 
convexity properties used in Eq. I|41ll do not hold in this 
case. Here thus we gives a lower bound for C\ obtained 
by assuming the encoding Eq. (|4*3^l . 



Ci ^ max \ (1 - 773P) 
pe[o,i] L 



Ho 



1 - (1 - m)p 



1 - mp 




}• (55) 



FIG. 3: Plots of the capacities (in bits per channel uses) of 
the channel T V1 , V2 of Eq. I|15p as a function of 772 € [0, 1 — 771] 
obtained for different 7/1 (this parameter increases from to 
1 moving along the arrows). In both graphics the dashed 
curve is the capacity of an amplitude damping channel with 
efficiency 1 — 772 which according to Eq. 14711 provides an up- 
per bound for the corresponding capacity of % )1 , rt2 - Part a): 
quantum capacity Q obtained by solving numerically the max- 
imization over p of Eq. 15311 . Notice that for 772 > 771 Q nullifies 
[35|. Part b), entanglement assisted classical capacity Ce of 

Eq. cm. 



IV. CONCLUSIONS 

In the previous section we calculated the capacities Q, 
Ce and C\ for the qubit amplitude damping channel T> v 
and for the channel 7~ Vl _ V2 . As discussed in Sec. Ill Al by 
identifying the parameters 77, 771 and 772 with the quanti- 
ties of Eqs. HA4(1 and (|A5|) we can use these results 
to analyze the correlations between distant points of the 
chain. A detailed analysis of the fidelity for state trans- 
fer J^Jl in the case of k = 1 and uniform coupling Jjy is 
given in Ref. . Such a paper solves the dependence of 77 
from the evolution time and from the length of the chain. 
In particular in 1] it is shown that the maximum value 
of 77 drops with the distance between the encoding spin 
A (at the site r) and the decoding spin B (at site s) as 
\r — s|~ 2 / 3 . Since according to Sec lIII Al the quantum ca- 
pacity Q vanishes when the transmissivity is ^ 0.5, this 
implies that long chains are not suitable to directly trans- 
mit quantum information using the scheme of Sec. Ill A II 
with k = 1. Even though the communication scenarios 
described by means of Eq. 0]) are incomplete (see the dis- 
cussion of Sec.[njl it is interesting to see in which way the 
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above difficulties can be overcome. One possibility is to 
optimize the values of rj by tailoring the interaction Jj j 
between the spins as proposed in Ref. 0. Alternatively 
one can use registers with k > 1 and then optimizing the 
value of rj of Eq. by means an appropriate choice Q 
of the coefficients Cj of Eq. (JSJ) . 

Yet a different approach would require the partition 
of the chain in smaller segments so to transfer informa- 
tion faithfully along the chain through swaps between 
neighboring segments analo gou sly to the ideas of quan- 
tum repeaters introduced in 36] . Breaking the chain can 
be achieved by applying locally external magnetic fields. 
A combination of a time dependent control of part of 
the chain together with perfect transmission of small seg- 
ments may lead to an improvement of the performances 
of spin chains to transport quantum information. On 
the other hand, the quantum capacity of the chain can 
be boosted by giving access the sender and receiver to a 
free two-way classical communication line |l8u37j . In this 
case, the ability of the channel in transmitting quantum 
signal can be increased by means of entanglement distilla- 
tion protocols and teleportation. For instance (as noticed 
also in [p) since any two qubits entanglement state is dis- 
tillable [33 , the two-way quantum capacity Q%.^ 7 of the 
chain is strictly greater than zero also when Q is null. To 
compute the exact value of Q2-wa y one needs to find the 
optimal distillation protocol for M. : unfortunately this is 
quite a challenging task. 



Acknowledgments 

This work was supported by the European Community 
under contract IST-SQUIBIT, IST-SQUBIT2, and RTN- 
Nanoscale Dynamics. 



APPENDIX A: DERIVATION OF EQ. itTH) 

The Hamiltonian (J2J preserves the total spin compo- 
nent along the z axis. The set of two-spin up states \j, I) 
is hence transformed into itself with unitary matrix ele- 
ments 

f^ rl ,{t)^{j,i\e- lHt/h W,l') ■ (Al) 

Using this transformation it is possible to show that at 
time t the spin chain state becomes 

N 

|*(t)>=aH>+/3 Y, AiMM > (A2) 

j>l=l 

where, for j > £ = 1, ■ ■ ■ , N, 

k 

AjA t )= d 3't'fjt,j't'(t) (A3) 

j'>i'=i 



are the time evolved of the coefficient dj t t of Eq. 
Equation (|A2() shows that only the fraction 

N 

m= £ \d,At)\ 2 (A4) 

j>l=N-k+l 

of \^(t)) has two spins up in B. Analogously 
N-k 

j>e=i 

is the probability of |^(t)) having both two spins up out- 
side of B, while 

N-k N 

m = i- m - m = £ E l d ^(*)l a . (A6) 

1=1 j=N-k+l 

is the probability of having one spin up in B and the 
other outside of B. By taking the partial trace of (|A2() 
over the first N — k spins of the chain we obtain the 
state of the quantum memory B of Eq. I|14|) . In such an 
expression, 

N 

\ti)B^ d jAtMt)/V*h- (A7) 

j>t=N-k+l 

is a rotation of \4>2}a that can be compensated for at the 
decoding stage. On the other hand, <t_b of Eq. (fHjl is a 
density matrix of B formed by states with one-spin-up 
vectors, i.e. 

N-k 

u B = \MQ) (MQ l/is. (A8) 

i=\ 

with 

JV 

\4>i{t)) B = Yl d u {m^)/VW) 

j=N-k+l 
N 

m (i) = Y \ d M\ 2 ■ (A9) 

j=N-k+l 

Notice that the rank of <jb is at most equal to the maxi- 
mum number of orthogonal one-spin up states of B, that 
is k. Moreover the support of <jb is clearly orthogonal to 
| JJ-) S and \<f)' 2 ) B - 

APPENDIX B: RELATION WITH THE BOSONIC 
LOSSY CHANNEL 

In the lossy Bosonic channel £ n [53, IH, |2^| an input 
Bosonic mode described by the annihilation operator a 
interacts, through a beam splitter of transmissivity 77, 
with the vacuum state \&)b of an external Bosonic mode 
described by the annihilation operator b. Any input state 
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p of the mode a is hence transformed by this map accord- 
ing to the equation 

£„(p) = Tt 6 [[/ (p®|0) 6 <0|)[/t] (Bl) 

where the trace is performed over the external mode b and 
where U is the beam splitter unitary operator defined by 

U^aU = y/Tj a+ y/l -r/b (B2) 
U f bU = y/rjb- y/l-rja. (B3) 

By restricting the inputs p to the Hilbert space spanned 
by the vacuum state and the one photon Fock state 
the map l|23[) has the same Kraus decomposition i|12|) 
of the qubit map T>, r Some capacities of the channel £ n 
have been solved under constrained average input pho- 
ton number: the classical capacity C is given in Ref. 
while the entanglement assisted capacity Ce and a lower 
bound for Q which is supposed to be thigh are given in 
Ref. [39j. Unfortunately, since an average input photon 
number constraint can not prevent the average message 
of £ v from being supported on Fock states with more 
than one photon (apart from the trivial case of zero av- 
erage photon number), the results obtained in [28l l39| 
provide only trivial upper bounds for the correspond- 
ing capacities of D^. For instance consider the entan- 
glement assisted case where the capacities of both the 
channels can be computed. In the input Hilbert space 
we are considering here, the average photon number of 
the transmitted message is provided by the average pop- 
ulation associated with the one photon Fock state. A 
fair comparison between the capacities of T)^ and £ v 
can be hence obtained by taking the value of Ce asso- 
ciated with a lossy Bosonic channel where the aver- 
age input photon number is given by the population p 
of Eq. I|38fl which maximizes the entanglement-assisted 
capacity of T> v . According to (3^ the capacity of £ v is 
then given by Ce = g{p) + g(vp) ~ <?((! — v)p)i where 
g(x) = (x + 1) log 2 {x + 1) — x log 2 x. A simple numerical 
analysis can be used to verify that this quantity is always 
bigger than the corresponding value (|38|l of T) v . 

The channels T)^ and £ v share many common features. 
In particular obeys to the same composition rule of 
given in Eq. 12 21) and it is degradable, since for any 
input p one has 

Sr,{p) = PS £i- v (fi) SP , (B4) 

where now P = e i7rb+b and S is the swap operator which 
transforms a in & and vice versa. This relation was used 
in |29j | without explicitly proving it, and applies to all 
Gaussian channels of the form l123[l where the external 
Bosonic mode b is prepared in a circularly symmetric 
input. For the sake of completeness, here we give an 
explicit derivation of Eq. (|B4|I in the case of the purely 
lossy Bosonic channel. 

Consider a generic input state p of the Bosonic channel, 
with characteristic function T(p) = Ti a [pD a (p)], i.e. 

p=J f^T(p)D a (-p), (B5) 



where D a {p) = exp^W — u*a] is the displacement op- 
erator of the input mode a, As shown in |2£j, the 
channel £ n transforms p into 

= / ^ F » > (B6) 

with 

r'(/i) = r(07/i)e" (1 ^ )lAl|2/2 - (B7) 

Analogously it is possible to verify that the complemen- 
tary map £ v (defined as in Eq. H25fl by replacing the 
partial trace over b with the partial trace over a in the 
Eq. IjBljl ) produces the following transformation, 

4(P) = / ^ f 'M D <>(-ri > (B8) 

where Db(p) is the displacement operator of the mode b 
and where 

f'(p) EE T(-^T^ p) e-"H 2 /2 . (B9) 

Suppose now r\ 0.5 and apply the lossy map £i- v / v 
to the state of Eq. (|B6|I : according to the composition 
rule (|22f) this will transform its symmetric characteristic 
function to 

r» ee r(yr^ p) e-^' 2 = f'(- M ) , (bio) 

producing the state 

s^(£, 1 ( P ))= / — r'(p)D a (p). (Bii) 
" J 7r 

This proves the identity (|B4|> since under the unitary 
transformation PS the annihilation operator a is trans- 
formed into — b and the state (|B11|I becomes equal to the 
output (|B8(1 of the composite map. 



APPENDIX C: SOME USEFUL RELATIONS 

In this appendix we provides some relations used in 
Sec. 11111 to derive the capacities of the qubit amplitude 
damping channel. 



1. Entropy of exchange 

The output state (V^ <g> t anc )($) associated with the 
purification |$) of Eq. (|2*§1) can be expressed in the com- 
putational basis { 1 00), 1 01), 1 10), |11)} of H A <8>H a nc (here 
|0) anc ee \Rq) while |l)onc is the component of |i?i) which 
is orthogonal to \Ro))- This gives the following 4x4 ma- 
trix 
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(1-^)7*^^ 



V^7 








|7| 2 
' 1— p 



'/ / 1-p 



»7 y 



(Cl) 



The matrix i|Cl(l has eigenvalues (two times degenerate) 
and 



A± = (l ± V(l - 2(1 - ^7)^)2 + 4(1 - 77)lTl 2 ) /2 
= A±(l-»7), (C2) 

with A± the eigenvalues ofT> v (p) given in Eq. (|32|l . The 
exchange entropy associated with p is hence given by 



S((V V <g> l a „ c )($)) 



-A + In A + - A_ In A_ 



= ff 2 (A+) , (C3) 

where, for cc € [0, 1], 

H 2 {x) = -x\og 2 x- (l-aj)log a (l-x) (C4) 

is the binary entropy Some useful relations of ^(z) 
are the following: 

1) The function Hi{z) is decreasing with respect to 
the variable 1 1/2 + z|. This property is a conse- 
quence of the fact that the entropy associated with 
a binary stochastic variable is maximum when the 
probabilities associated with different outcomes are 
equal. 



2) The function H 2 ((l + \/l - z 2 )/2) is convex with 
respect to z. This property can be easily verified 



and is related with the convexity of the entangle- 
ment of formation with respect to concurrence in 
qubit systems |42^ . 



2. Dependence on 7 of the coherent information 

As discussed in the text, the quantum capacity Q of 
T> v is obtained by maximizing the function J(p, I7I 2 ) of 
Eq. over all values of p G [0, 1] and I7I 2 G [0, (l-p)p\. 
For 77 1/2, this expression is decreasing in I7I 2 , i.e. 
for any \{3\ G [0, 1] it achieves the maximum value for 
7 = 0. The proof of this result is quite tedious, but can 
be obtained analytically by studying the partial deriva- 
tive of with respect to the parameter |7| 2 . We skip 
all the details of this analysis which is not of a fundamen- 
tal interest, and simply observe that the problem can be 
reduced to studying the properties of the function 



1 - - 2x) 2 +Axy 
1 + - 2.t) 2 + 4xy 



x/y/(l-2x) 2 +4xy 



(C5) 



on the domain x G [0, 1 — y], for any y G [0,1]. One 
can then verify that f y (x) is decreasing in x: this guar- 
antees that J(p, I7I 2 ) is monotonically decreasing in |7| 2 , 
yielding the thesis. 
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